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We study how to control the dynamics of tripartite entanglement among optical cavities using 
non-Markovian baths. In particular, we demonstrate how the reservoir engineering through the 
utilization of non-Markovian baths with different types of Lorentzian and ohmic spectral densities 
can lead to an entanglement survival for longer times and in some cases considerable regain of 
seemingly lost entanglement. Both of these behaviors indicate a better sustainability of entanglement 
(in time) compared to the usual Markovian bath situations which assumes a flat spectrum of the 
bath around the system resonant frequency. Our scheme shows these effects in the context of optical 
cavities starting off in a maximally entangled W and Greenberger-Horne-Zeilinger (GHZ) tripartite 
states. In Lorentzian cases we find that the far detuned double Lorentzian baths with small coupling 
strengths and for ohmic type baths super-ohmic environments with smaller cutoff frequencies are 
the best candidates for preserving entanglement among cavities for significant amount of time. A 
non-Markovian quantum jump approach is employed to understand the entanglement dynamics in 
these cases, especially to recognize the collapse and revival of the entanglement in both W and GHZ 
states. 


I. INTRODUCTION 

The problem of maintaining entanglement SI] among 
quantum systems that are open to their environments 
is one of the most challenging issues in quantum in¬ 
formation processing today [2|. It is well known that, 
even if the interactions are weak, the environment 
(bath/reservoir) dramatically alters the dynamics of the 
quantum systems, which results in an irreversible loss 
of information from the system to the environment: a 
phenomena more commonly known as the decoherence 
[3|. As a result, this lost information doesn’t play any 
further role in the dynamics of the system and the 
system loses its pure quantum features (in particular 
entanglement). 

The successful generation and longer survival of en¬ 
tanglement is the central ingredient to a variety of 
quantum information processing tasks such as quantum 
teleportation [3], quantum dense coding [5], quantum 
cryptography [6] and quantum computing [7j. Entan¬ 
glement also plays a crucial role in the foundations of 
quantum mechanics and in the non-classical behaviour 
of correlated many-body quantum systems [5]. Based 
on this broad range of applications, devising novel 
ways of preserving entanglement among open quantum 
systems (for considerable times) is of central importance 
for the practical applications of multipartite quantum 
systems in future quantum technologies. In order to 
control and reduce the effects of decoherence, in recent 
years reservoir engineering techniques are developed 
at a practical level, which are now revising the role 
of the system-environment interactions PHI]. These 
techniques are leading to new ways to manage the 
properties of environment which can help in sustaining 
the pure quantum features of the system for extended 
times. 


Motivated by these techniques, in this work we ad¬ 
dress this decisive question that once the entanglement 
among open quantum systems is established, how we 
can control/save it for a sufficient amount of time so 
that the quantum protocol (related to the quantum 
communication or quantum computation under con¬ 
sideration) can be successfully performed. We answer 
this question in the context of tripartite entanglement 
which is assumed to be generated among three open 
and independent optical cavities. Our choice of studying 
tripartite entanglement is based on the fact that, 
despite of intense research in entanglement for last two 
decades, at present our full understanding regarding 
entanglement is merely limited to bipartite systems. 
Although there are many criteria introduced to quantify 
multipartite entanglement of mixed and pure states 
jT], a full framework is still lacking. One particular 
class of multipartite entangled states which is relatively 
easy to study is the tripartite entangled states (12]. W 
and GHZ states are the celebrated examples belonging 
to this class of entangled states. Along with many 
applications in quantum information processing (as in¬ 
dicated in the references above) the tripartite entangled 
states may also provide a platform to understand the 
differences and similarities in the behavior of multipar¬ 
tite and bipartite entanglement at the fundamental level. 

The basic idea that we are going to exert in the present 
study will be to employ the technique of reservoir en¬ 
gineering through the utilization of non-Markovian 
baths (with different spectral profiles) for freezing the 
entanglement among tripartite-cavity system. More 
precisely, we’ll investigate the type of non-Markovian 
baths that are most suitable for the entanglement 
preservation from environmental decoherence. In recent 
years non-Markovian baths have attracted considerable 
attention due to their diverse applications in the foun- 
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dations of quantum theory, in quantum computation 
and quantum communications m I14j . The reason for 
utilizing non-Markovian baths for entanglement survival 
purposes, lies in the fact that in the standard approach 
of studying open quantum systems, environments are 
usually assumed to be memory-less (Markovian) on 
the time scales of the system dynamics. As a result, 
environment always act as a sink for the information 
contained in the system. It is well known that the 
bipartite entangled systems suffer entanglement sudden 
death (ESD) dSUISj due to this effect. 

Such a fast decay of entanglement is not desirable for 
several quantum information and quantum computation 
protocols where entanglement plays a key role m- 
Hence we employ non-Markovian environments where 
we’ll have more freedom of choosing environments 
with different spectral functions as well as we can 
identify different ways of probing the environment and 
environment-system interaction such that entanglement 
can be sustained for longer times. In here, we are 
going to consider the non-Markovian environmental 
models which have experimental relevance too, such as 
in condensed matter [15] . biological m and quantum 
optical [201 HU systems (to name a few). Except one 
case of sub-ohmic spectral function, we find that using 
non-Markovian baths with different Lorentzian and 
ohmic type spectra can better perform in the longer 
entanglement survival compared to usual Markovian 
baths. In particular, weakly interacting far detunned 
double Lorentzian baths and super-ohmic baths with 
smaller cut-off frequencies establish an entanglement 
which is more robust against environmental losses. 
Moreover, a non-Markovian quantum jump approach 
(NMQJA as introduced by Piilo et.al [22]) is found to be 
sufficient to explain the loss and gain of entanglement in 
terms of positive and negative cycles of time-dependent 
decay rates. Here we would to mention that the bipartite 
and tripartite qubit systems dynamics when coupled 
to Markovian and non-Markovian baths is already a 
studied subject [2M251 . The main novelty of our work is 
the consideration of various environmental models (with 
finite temperratures) in this context and the application 
of NMQJA to understand the entanglement dynamics, 
which to our knowledge has not been reported previously 
for the system under consideration. 

The paper is organized as follows: in Sec. II, we model 
the system by writing down the Hamiltonian and include 
the dissipative dynamics of the system by coupling the 
cavities to their respective baths with different spectral 
densities. In Sec. Ill, we present results and discus¬ 
sion where we’ll first consider the entanglement dynamics 
when the optical cavities (initially either in W or in GHZ 
state) are coupled to a Markovian bath. Effect of photon 
hopping between cavities will also be considered. This 
section will be useful for later discussions on a compari¬ 
son between entanglement dynamics for Markovian and 


non-Markovian baths cases. Next we study the effect of 
non-Markovian baths on the cavity based tripartite en¬ 
tanglement dynamics of W and GHZ states. Finally in 
Sec. V we’ll summarize the results and report conclusions 
of this work. 


II. HAMILTONIAN AND DISSIPATIVE 
DYNAMICS 

System under study comprises of three independent 
and open optical cavities. For simplicity it is assumed 
that each cavity has a single isolated resonant mode, 
which are given by frequencies w c i,w C 2 and uj c 3 in the 
first, second and third cavity respectively. Destruction of 
photon in these cavities is described by the annihilation 
operator cq Vi = 1,2,3 (throughout the article we’ll use 
same range of values for index i). Under rotating wave 
approximation, the total/global Hamiltonian describing 
the system Hamiltonian (Hg ), bath Hamiltonian ( Hb) 
and system-bath interaction Hamiltonian ( Hsb ) is given 
by: 


H sys — Hg + Hb + Hsb — ^ ^ b ki b ki 

i= 1 V ki 

+ '^2,hg ki {a i b\ i +a\b ki } 

ki J 

(1) 

Baths are assumed to have a continuum of modes while 
b ki operator annihilates a photon in the /cth mode 
(frequency oj ki ) of the ith bath. The system (cavity) 
mode interaction with the bath modes is characterized 
by the real coupling rate g ki without loss of generality. 
Non-vanishing commutation relations are given by: 
[a,,a]] = Sij , [b ki ,bl.] = = 1,2,3 and j = 1,2,3. 

When the cavities are coupled with non-Markovian 
baths, the system dynamics can be described in terms of 
the so-called Exact Master Equations (EME) [22]. Basic 
assumption in these types of Master equations is the 
dependence of the time derivative of the system density 
operator on the state of the system at the present time. 
Although the general structure of EME resembles the 
Lindblad Master equations but due to the presence 
of time-dependent functions (which are related to the 
photon decay rates and energy shifts m\) the EME are 
generally non-Lindblad per se. It turns out that these 
time-dependent decay rates are responsible for the pro¬ 
duction of memory effects (non-Markovianity) in EME 
(we’ll notice this fact in the context of entanglement 
evolution in later sections). 

In interaction picture, EME for present systems takes 
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the following form 


dp s {t) 


oti{t) a\p s {t)di - p s {t)did\ 


dt 

i—1 

+ a*(t)\ d\p s {t)ai -did\p s {t) ) + Pi(t)[ dip s {t)d\- 


a\dip s it) ) + P*(t) ( dip s {t)d\ - p s {t)d\di 


( 2 ) 
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The time-dependent functions ai{t) and Pi{t) are re¬ 
lated to the non-Markovian baths’ spectral densities 
Ji(oj) in the following way: 

pt pOO 

Oi{t) = dt i / dwJ l (w)]V(a;)e l( ^-“ )(t -* l) (3a) 

Jo Jo 

pt poo 

/3i(t) = / dt ± / duJi{u)(N(u) + (3b) 

do do 

while bath spectral densities Ji{ui) are in turn dependent 
on the cavity-reservoir interaction rate through 
Ji(w) = f dt Y^ k J3D]. For the sake of 

simplicity all baths are assumed to have same average 
number of photons -/V(w) as given by the Planck’s dis¬ 
tribution ehul/ k B T —i while ks is the Boltzmann constant 
and T is the absolute temperature of the bath. 

Choosing the basis set according to the truncated 
Hilbert Space (H) relevant to the present three cavity 
system: {|000> ->• |1), |100) -> |2), |010) -> |3), |001) -> 

14), |HO) |5), 1101) -> 16), |011) 17), |111) -»■ |8)} 

(while the first, second and third slots in the ket 
represent the number of photons in the first, sec¬ 
ond and third cavities respectively), and assum¬ 
ing all the time-dependent functions to be the 
same i.e. a\(t) = ct 2 {t) = 0:3 (t) = a(t) and 

Pi (t) = @2 {t) = /?3 (t) = P{t), one can calculate the 
time evolution of density matrix elements. For W state, 
initially non-vanishing density matrix elements are: 
P22it = 0) = P23it = 0) = P24)it = 0) = P 32 (t = 0) = 
P33it = 0) = p 34 (t = 0) = p i2 it = 0) = P43it = 0) = 
P 44 ^t = 0) = 1/3. For GHZ state initial conditions are: 
Pn{t = 0) = pisit = 0) = psiit = 0) = P88it = 0) = 1/2. 
Full solution of the density operator p s (t) is expressed 
in terms of following matrix: 


CJ 

FIG. 1: Spectral density functions describing 
different environmental models. Parametere used 
are: For Loretzian functions (single, double and 
band-gap Lorentzian) a 4 = a^i = oll 2 — 2, 
r = r 4 = r 2 = 1, 0 J bc = 20, W D1 = W D2 = 1/2 and 
Wbi = 2, W B 2 = 1- For ohmic type densities we 
have used u cut = 15,10, 2 and a = 1,0.6, 0.1 for 
sub-ohmic, ohmic and super-ohmic functions 
respectively. Note that we have also plotted the flat 
spectrum density which is commonly assumed in 
the case of Markovian baths. 


(and experimental relevant) types of environmental mod¬ 
els having: single Lorentzian (Jsl), double Lorentzian 
(Jdl), band-gap Lorentzian (Jbl), sub-ohmic (Jsbo), 
ohmic (Jo) and super-ohmic (Jspo) spectral densities 
as expressed below and shown in Fig. 1: 

JsL<u> = S (^-^J + (r ftr (5a) 

Jdl( ") = Wdi (irF^yVow) 

+WD1 ( y (5b) 
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In general, non-Markovian baths can have any type of 
spectral profile. In present study, we’ll consider different 


J( W ) = a Wc 1 -/ W s e-/^. (5d) 

The non-Markovian bath with a Lorentzian spectral 
density can be experimentally realized by coupling the 
quantum system with a high finesse (leaky) optical cav¬ 
ity (which we’ll call the bath cavity). In pseudomode 
picture of non-Markovian dynamics [31] system’s cav¬ 
ity mode couples with the pseudomode of the bath cav¬ 
ity which in turn couples with the outside Markovian 
environment. In single Lorentzian density, c*i, T and 
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uibc express the coupling between the bath cavity and 
the outside environment, width of the spectral distri¬ 
bution and the resonant frequency of bath cavity re¬ 
spectively. In double (band-gap) Lorentzian we have 
summed (differed) two Lorentzians with the parame¬ 
ters (describing each Lorentzian function) having the 
same meaning as in the case of single Lorentzian model. 
In addition to these parameters, in double Lorentzian 
model both Lorentzians are weighted by the condition: 
Wm + Wd 2 = 1 and in band-gap Lorentzian model the 
conditions Wbi — Wb 2 = 1 and T 2 < lb are respected 
to ensure positivity of the spectral density as a func¬ 
tion of u>. To examine a bath which considers an even 
broader range of frequencies, in last example (Eq.5(d)) 
we have summarized sub-ohmic, ohmic and super-ohmic 
spectral densities, by introducing spectral parameter s. 
Note that in this equation by choosing s = 1/2,1 and 
s = 3, we can obtain J(w) —> Jsbo(w), J(w) —> Jo(w) 
and J{oj) — Jspo(w) respectively. To cut off the higher 
frequencies (to avoid divergences) u> cu t is introduced in 
the spectral density J(oj) while a represents a dimension¬ 
less coupling constant. 


III. RESULTS AND DISCUSSION 
A. Entanglement dynamics under Markovian baths 

1. Without photon hoping 


Note that one can also obtain the above Master equa¬ 
tion directly from Eq.[2] by setting Ji(co) = Hi/2^ and 
Ni(u) = Ni = e no Jci/ l B T_ 1 - The decay rate of energy 

from the ith cavity is denoted by Ki and N identifies the 
temperature of the heat bath which is assumed to be 
same for all baths. The density matrix p(t) (describing 
the state of the optical cavities) can be used to quan¬ 
tify the amount and temporal evolution of entanglement. 
There are several different measures of entanglement that 
can be used to quantify bipartite and multipartite entan¬ 
glement pQ. One measure which we are going to use 
here is called the negativity (AT), which follows from the 
Peres-Horodecki Positive Partial Transpose (PPT) crite¬ 
rion and it can be calculated for systems of arbitrary 
Hilbert space dimensions. The range of A f lies between 
0 and 1 where 1 refers to maximum entanglement and 0 
to completely separable states. Negativity is defined as: 

A f{t) = max ^0, —2 E A *) ( 7 ) 

where the sum is taken over the negative eigenvalues \ 
of the partially transposed density matrix pf T . Partial 
transposition is taken with respect to one of the cavi¬ 
ties only and here we’ll perform it with respect to the 
first cavity. We find that (independent of bath temper¬ 
ature ( N )), among eight eigenvalues of pf T always only 
one eigenvalue turns out to be negative. For an initial 
preparation of the system in W state: 


Different experimental techniques and theoretical pro¬ 
posals are developed for the successful generation of tri¬ 
partite entanglement among quantum systems. Four 
wave mixing in cold atomic gases [32j , spatial mode para¬ 
metric down conversion [33], anisotropic exchange inter¬ 
actions [MJ, single spins in diamond )3!> and scheme 
based on cavity quantum electrodynamics (CQED) (pho¬ 
tons passing through optical cavities) [36., are few promi¬ 
nent examples. 

To highlight the differences between Markovian and non- 
Markovian entanglement dynamics, in this section we ad¬ 
dress the Markovian bath case when three cavity system 
is initially prepared in a tripartite entangled state (either 
W or GHZ state). For simplicity we’ll begin by assuming 
no direct or indirect coupling among cavities and among 
baths as well. With a single photon restriction in each 
cavity and following the standard reservoir theory [25] . 
the dynamics of such a two-level three-cavity system can 
be described by the well known Lindblad Master equa¬ 
tion: 

df) ^ = E (' K i( N + ^){aip s {t)a | - | a\aip s {t)- 

2 = 1 \ 

^p s (t)a\ai) + KiN(a\p s {t)ai - ^aidjp s (t ) - ^p s (t)dia\) 

( 6 )' 


mt = 0)} w = Eflioo) +1010) +1001)) (8) 

we find the following analytic form of Negativity: 


A/w = max 


0,(e- Kt -l)(l + (l + ^e- 2Kt (e" Kt 
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And for a GHZ initial state: 

l^ = 0)} G *z = ^(|000> + |lll» 

Negativity turns out to be: 


I)” 2 )) 

(9) 

( 10 ) 


Nghz = max 


0, de~ Kt (e~ Kt - 1) + i e - 3Kt (4e 3Kt + 

o o 


-(4 - l2e Kt + 13e 2Kt - 10e 3Kt + e 4Kt ))) 

O 


( 11 ) 


In FIG. 2(a),-(b) we have plotted Negativity for both 
W and GHZ states with varying values of N. In both 
cases we notice an asymptotic and irretrievable loss of 
entanglement due to the cavities interaction with the 
Markovian baths. This decaying behaviour becomes even 
steeper with increasing the bath temperature ( N ). In the 
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(a) (b) (c) 



FIG. 2: Time evolution of tripartite entanglemnet among optical cavities characterized by negativity, when 
cavities were initially prepared in a maximally entangled (a) W and (b) GHZ state excluding the effect of 
photon hoping, (c) W-state entanglement evolution at zero temperature when a two-way strong photon hoping 
among cavities is included (£12 = £23 = 5k with k\ = k 2 = K3 = k) . Note that we have shown the result only for 
the case of W-state, but the main points of this plot appear to be the identical for an initial GHZ state as well. 


case of two qubits (bipartite systems) coupled to their 
independent reservoirs such a decaying profile is a well 
known phenomena and it is regarded as the ESD m- 
We notice that a similar decaying behaviour of entangle¬ 
ment extends down to tripartite entangled cavities also 
in the presence of Markovian environments. 

2. Effect of Photon Hoping 

In order to retain entanglement for longer times in the 
presence of usual Markovian baths, one proposal which 
can be adopted is to include the effect of photon hop¬ 
ing among cavities. Such processes can be made more 
probable by placing cavities closer so that the photon 
wavepacket can overlap/tunnel to the other cavity. Al¬ 
though such a system will not be very useful for the pur¬ 
poses of long distance quantum communications (such as 
in continuous variable based quantum key distribution 
protocols [38] ), but still it may be used to store entangle¬ 
ment for extended times. With the inclusion of photon 
hoping, the system Hamiltonian as given in Eq.(Jl]) ex¬ 
tends to: 

Hnew = H sys + h£ 12 (a\a 2 + a\a 1) + hf, 23 (a\a 3 + d^a 2 ) 

( 12 ) 

where H sys is given in Eq.(l) and the rest of terms 
on right side of above equation are hoping terms (£12 
and £23 describing the hoping strength between cavity 
1 and 2 and cavity 2 and 3 respectively.) In FIG.(2)c 
we have plotted the entanglement evolution for W-state. 
We notice that due to the hoping of photon, J\T(t) shows 
slight oscillations while decaying. Along with this it 
appears that the negativity vanishes at the same time 
(around t = 2k -1 ) as it vanished for no hoping case (see 
FIG. (2)a). But a close inspection reveals the presence 
of very small entanglement for later times t > 2 k _1 (as 
shown in the inset of the figure). This small entangle¬ 
ment exhibits the possibility of photons going back and 


forth between cavities due to strong hoping before being 
leaked out. 

Although such a proposal can be used to save the entan¬ 
glement in the system for longer times but the amount 
of entanglement left in the system turns out to be so 
small that it can hardly be useful for any practical ap¬ 
plication. Keeping in view this fact, in next section we’ll 
investigate the effect of non-Markovian environments on 
the entanglement survival among cavities. 


B. Non-Markovian baths 

There are many quantum information protocols 
that however, require the sustainibility of considerable 
amount of entanglement for longer times. Clearly cou¬ 
pling of quantum systems with Markovian baths (even 
with the possibility of photon hoping) will not support 
these protocols and hence novel approaches are needed 
for the entanglement storage for desired amount of time. 
In order to address this problem, in recent years the 
method of entanglement control by using quantum feed¬ 
back networks is developed j35]. We on the contrary 
here propose to exploit the the technique of reservoir en¬ 
gineering through the utilization of non-Markovian baths 
for this purpose. 

We discuss the dynamical behavior of tripartite entangle¬ 
ment, when cavities are coupled to non-Markovian baths 
of different Lorentzian and ohmic type spectral profiles. 
Here we’ll first consider a relatively simple situation when 
the temperature of all baths is set to zero. Later on (in 
section 3.3) we’ll also consider the effect of finite temper¬ 
atures. The main aim here will not only to compare the 
effect of Markovianity and non-Markovianity on entan¬ 
glement but also to analyze which type of spectral den¬ 
sity of the non-Markovian baths can surpass other types 
of densities in producing entanglement which is most ro¬ 
bust against environmental decoherence. 
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1. Lorentzain type Profiles 

We’ll start with the examples of non-Markovian baths 
with Lorentzian type spectral densities. For baths at zero 
temperatures, EME takes a simpler form with afit) = 
0 and 0i(t) = J* dh The 

time-dependent decay rate «j(i) is related to the time- 
dependent function p j(f) through the relation «,(£) = 
2Re[/3;(t)]. For the example of single Lorentzian bath 
(characterized by density Jsl{ w)), time-dependent de¬ 
cay rate takes the form: 


(13) 

whereas for double and band-gap Lorentzians we’ll con¬ 
sider two single Lorentzian decay rates nf L1 (t) and 
k f L 2 (t) with different widths of distribution (Fi and T 2 ) 
which are weighted appropriately and added/subtracted 
to give the following time-dependent decay rates: 

' L {t) = W D1 Kf L1 (t) + W D2 Kf L2 (t) (14a) 

nf L (t) = W B1 K? L1 (t) - W B2 nf L2 (t) (14b) 

Si = (i u> c i — cube) is the detunning frequency and for sim¬ 
plicity we have assumed that all time-dependent decay 
rates are the same i.e. Ki(t) = n 2 (t) = K 3 (t) = K(t). 

• On resonance case: With these time-dependent 
functions, we numerically solve the equations of motion 
of the density matrix elements (as outlined in Appendix 
A) and then the Af(t) is calculated by applying the 
PPT criterion on p s (t). Similar to Markovian bath 
case, we notice that only one of the eigenvalues in the 
partial transposed density matrix (we call it A_) exhibit 
negative character as required by the definition of 
negativity, thus leading to: Af(t) = max( 0, 2A_(f)). In 
Fig. 3(a),-(c) we have plotted the dynamics of tripartite 
entanglement when cavities are initially prepared in W 
and GHZ states respectively. All three examples (single, 
-double and -band-gap Lorentzians) are shown along 
with Markovian case. Here firstly we have presented the 
results for on resonance case i.e. 5 = 0. 

In Fig. 3(a),-(c), we notice that utilizing a non- 
Markovian bath of any type of Lorentzian spectral profile 
results in somewhat slower decay (and hence longer sur¬ 
vival) of entanglement compared to the usual Markovian 
bath case. Double Lorentzian bath makes this behav¬ 
ior most pronounced. GHZ state has initially (at t = 0) 
slightly more entanglement than W state and that dif¬ 
ference prevails throughout the plot. In W and GHZ 
states, system start off in maximally entangled tripartite 
states, but reaching at £ ~ 8,5,4 and 3 (units of time) for 
double, -single, -band-gap Lorentzian and flat Markovian 
baths respectively, entanglement completely dies out due 
to photon leakage. 


The different decaying slopes of A f(t) in all curves, can be 
explained by plotting the corresponding time-dependent 
decay rate (as shown in Fig. 3(b)). Since the decay 
rates have time-dependence therefore we calculate the 
average decay rate over entire time interval (i.e. «(£) = 

St ft +AT whereas in present case AT = 10 w” 1 ). 

We find that for double -single and band-gap Lorentzians 
the average decay rate turns out to be 0.42, 0.74 and 1.37 
wfi 1 respectively. These time-dependent rates represent 
effective decay seen by the cavities and it explains why 
for example the entanglement decays most slowly (fastly) 
in double (band-gap) Lorentzian case compared to other 
situations. 

• Off resonance case: In Fig. 4 we plot the negativ¬ 
ity and time-dependent decay rates as a function of time 
under finite detunning case. We note that by introducing 
5^0, entanglement start to show oscillations superim¬ 
posed by an overall decay. Why entanglement shows fast 
and slow decay and collapse and revival against different 
spectral baths for off resonance case? We’ll answer this 
question by employing the NMQJA developed by Piilo 
et.al in [22] (see Appendix B for a brief review). 
Explanation of negativity’s behavior based non- 
Markovian quantum jump approach (NMQJA) 
To apply NMQJA, we have plotted the time-dependent 
decay rates for different types of Lorentzian environmen¬ 
tal models in all cases. We note that the periodic tem¬ 
poral structure of decay rate translates to A f(t). During 
the positive cycle of decay rate, entanglement decreases 
(positive jump occurs), while during the negative cycle 
of the decay rate, non-Markovian environment transfers 
(some part of) entanglement back to the system (nega¬ 
tive jump occurs). 

It seems surprising though that the entanglement evolu¬ 
tion for both W and GHZ states behave in exact same 
manner. In terms of NMQJA, this difference can be seen 
by considering for example that the system start off in 
a GHZ state. Say during the positive cycle, one photon 
is leaked by each cavity then the initial state of cavities 
| 111) can be in any one of the kets |011) , 1101), 1110). But 
it is also possible that during this time another photon 
is emitted then the system will be found in either 1100) 
or in 1010), 1001). Now later when negative cycle arrives 
two negative jumps will be needed to obtain the initial 
maximally entangled GHZ state. Clearly for an initial 
W-state only one positive jump is possible which can be 
compensated by one negative jump. Thus one will expect 
less entanglement regain for GHZ state as compared to 
W state upon a negative cycle of decay rate, which is of- 
course not the case (as seen in Fig. 4). This identicalness 
of Af(t) behaviors in W and GHZ states can be explained 
by assuming that the decay rate is weak ({T!^} < a) 
and hence the possibility of a second photon decay is al¬ 
most negligible. At a single photon leakage level thus one 
will observe an identical Af(t) profile (as seen in Fig. 4). 
In single Lorentzian case (Fig. 4(a) and (d)) k(£) os¬ 
cillates between 0.58 and -0.28 and hence entanglement 


1 — e 
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(a) (b) 



(c) 



FIG. 3: Temporal profile of entanglement (negativity) among three optical cavities when system start off in a 
(a) W state (c) GHZ state, for different Lorentzian non-Markovian baths as well as for flat Markovian bath. In 
all plots 5 = u> c — ujbc = 0 and for single Lorentzian bath ol = 2, F = 0.1u; c . For double and band-gap 
Lorentzian common parameters are: chli = aL 2 = 2,Fi = 0.1w c ,r2 = 0.01w c whereas for double (band-gap) 
Lorentzian Hq = Wi = 1/2 (W\ = 2, W-j = 1) are chosen. 


shows clear collapses and revivals. While for double 
Lorentzian case (Fig. 4(b) and (e)) the amplitude of 
the decay decreases and remains between 0.3 and -0.2, 
thus Af(t) doesn’t fully dies out before the revival cy¬ 
cle arrives and entanglement again shows an increase. 
Finally we remark that the double Lorentzian bath pro¬ 
duces the most robust entanglement against decoherence 
(considerable entanglement survives till 20re -1 ) dominat¬ 
ing all other types of baths considered in Fig. 4. Here we 
would also like to point out that there has been some 
studies conducted in past where effect of detuning on the 
entanglement preservation between two qubits has been 
investigated [4(71 l4l] . It is interesting to note that the 
detuning turns out to be a useful parameter for entan¬ 
glement storage for the tripartite system under study as 
well. 

• Effect of decreasing coupling constants: In 

order to further elevate the tripartite entanglement for 
initial times, we varied different parameters involved in 
the dynamics of negativity. We identified that by de¬ 
creasing the strength of dimensionless coupling constants 
(a£, Oil, 0 ^ 2 ) this issue can be resolved. Fig. 5 exhibits 
this effect, where we have selected ol(= oli = &L 2 = 2), 
three times smaller then the values considered in Fig. 4. 
As a result, non-Markovian entanglement now shows 
much higher values than the one found in Fig. 4 for ini¬ 
tial times in particular and for entire time range in gen¬ 
eral. From the perspective of NMQJA this happens due 
to smaller amplitudes of the decay rate k(£) compared 
to the ones seen in Fig. 4. Alternately in the pseudo¬ 
mode picture of non-Markovian baths this behavior can 
be explained by arguing that with the consideration of 
smaller coupling constants the interaction of the bath 
cavity mode with the outside Markovian environment 
decreases as a result photon remains inside the system 
cavity and A f(t) shows longer survival. 

• Far detunned case: Another range of parameter 
that we found relevant is the regime where bath cavity 
is considered to be far detunned from the system cavity. 


We find that (as shown in Fig. 6) the consequences of 
considering far-off resonant Lorentzian baths (<5 = 5 u c ) 
is substantial. In all types of Lorentzian baths entan¬ 
glement shows tiny oscillations with no more than 20% 
of decay to the initial value upto t = 20 oj~ 1 . Especially 
for double Lorentzian case, Af(t) remains within 10% of 
its initial value, hence showing very small deterioration. 
From the perspective of NMQJA, this behavior is at¬ 
tributed to the smallest amplitude gained compared to 
cases considered in Fig. 4 and 5. Thus we concluded 
that the far detunned Lorentzian baths produces an al¬ 
most decoherence free evolution of entanglement among 
cavities irrespective of the initial entangled state. 


2. Ohmic type Profiles 


Next we address the non-Markovian baths with ohmic 
type of spectral densities. Such type of spectral den¬ 
sity corresponds to classical models of velocity-dependent 
friction forces and hence are termed as the ohmic func¬ 
tions. For the density function given in Eq.5(d), the time- 
dependent decay rate in present case turns out to be: 


, „ ot 
= - 


1 — cos[(s — l)tan 1 (u cu tt)\G(s — 1) 


(i 


t t 2 r 


(15) 


where G(s — 1) is the Euler gamma function. To in¬ 
vestigate the influence of these types of non-Markovian 
models on tripartite entanglement among cavities, in 
Fig. 7 we have plotted the negativity dynamics and time- 
dependent decay rates for sub-ohmic, ohmic and super- 
ohmic baths. For both (W and GHZ) type of initially 
entangled cavities, we find that due to smaller cut-off 
frequencies super-ohmic baths work very well compared 
to the Markovian and other types of ohmic baths to store 
entanglement (as shown in Fig.7(a),-(b) where consider¬ 
able amount of A f(t) survives until t ~ 10). Utilization of 
sub-ohmic and ohmic baths on the other hand turns out 
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(c) 
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FIG. 4: Entanglement and time-dependent decay rates plotted versus dimensionless time io c t for Lorentzian 
models for an initial (a),(b),(c) W and (d),(e),(f) GHZ state inclusing detunning S = lu c . Other parameters 
used (in all parts of the figure) are: a = a\ = ot\ = 6, T = Ti = 0.1w c ,r 2 = 0.01w c . Weightage factors (Wi, W 2 ) 
are the same as chosen in Fig. 3. Note that due to non-zero detunning, entanglement shows oscillatory 
character, with a possibility of collapse and revival. 




FIG. 5: Effect of reducing the dimensionless coupling constant cx.L,cx.Li,a-L 2 on non-Markovian tripartite 
entanglement dynamics. In all plots = «li = «L 2 = 2 is considered and time axis is measured in u>~ x units, 
rest of the parameters are the same as used in Fig. 4. Small amplitudes and oscillations in n(t) explains the 
higher and oscillatory negativity. 


to be not very useful. In ohmic case although the neg¬ 
ativity remains slightly higher than the Markovian case 
entanglement, but it dies out almost at the same time ( 
t ~ 2) and after that due to the absence of negative cycle 
in k (t) no revival occurs. Most strikingly, the sub-ohmic 
non-Markovian baths works even worse than the Marko¬ 
vian baths for the purposes of sustaining entanglement 
and Af(t) vanishes almost at t ~ 0.5. Again time average 
decay can give an effective rate seen by each cavity which 
is responsible for smaller/faster negativity decay. 


C. Finite temperatures 


The theoretical model described above can easily be 
extended to incorporate the consequences of finite tem¬ 
peratures. For simplicity here we’ll present our results for 
the case of non-zero temperature Lorentzian types of non- 
Markovian baths at low temperatures i.e. Huj c i > ksTi = 
1//3;. For the Lorentzian types of baths we find that 
under this temperature regime, time-dependent function 
(with assumption that a\(t) = a 2 (t) = (t) = a(t) and 
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FIG. 6: Tripartite-cavity entanglement variation in time for far detunned case 5 = 5 uj c (rest of the parameters 
are the same as used in Fig. 4 and with time measured in w” 1 units). Exploiting the off-resonance photonic 
decays, we notice that N(t) shows an extremely robust behavior against environmnetal interactions. Substantial 
entanglement survives upto t ~ 20co ~ 1 in all Lorentzian cases which is due to very small oscillation amplitude 
gained by n{t) (as shown in Fig. 6(b)) in accordance to NMQJA. 




FIG. 7: Plot of A f(t) versus t when an initial (a) W and (c) GHZ states are exposed to sub-ohmic (s = 1/2), 
ohmic (s —> 1) and super-ohmic (s = 3) non-Markovian baths. Note that at exactly s = 1 the time-dependent 
function n(t) blows up so to consider ohmic case we have taken the limit s —> 1. Other parameters used are: 
o-sub = 0.1 ,oj cu t = 2 for sub-ohmic, a„/, m = 0 ,6,uj cut = 10 for ohmic and a sup = 1 ,w cut = 15 for super-ohmic 
models and time axis is measured in a~^ p units. Figure inset shows time-dependent decay rate for super-ohmic 
bath. 


/3i(f) = /? 2 (t) = Ps(t) = /3(f) ) take the following form: 


a{t) = 


a L T 2 N 

2(<5 2 +r 2 ) 



a L T 2 Ne~ ihl3T 
2F(<5 + iT) 



P(t) = a{t) + 


a L V 

2 


1 - e~ rt e~ iSt \ 

r + iS J’ 


(16a) 

(16b) 


Corresponding time-dependent decay rate «(t) which 
for the simple case of zero detunnings (5 = 0) turns out 
to be: 

n(t) = ql( 1 — e~ rt )(lVcos(?i/3r) + 1) (17) 

With these time-dependent functions we plot the nega¬ 
tivity (for both initial W and GHZ states) in Fig. 8. In 


all plots we have taken N = 0.1. Compared to Fig. 3 
(in which N = 0) we find that main effect of non-zero 
temperature is an enhanced decays of in all single, 
-double and -band-gap Lorentzian models. For example, 
in double Lorentzian case AT(t) now dies out at around 
~ 5W” 1 while in N = 0 it died at ~ 6W” 1 . Rest of the 
features including the overall temporal profile, order of 
highest to lowest prolonging time periods of Af(t) in dif¬ 
ferent Lorentzian models and no-revival behavior remains 
the same as noticed in IV = 0 case. Most interestingly, 
when N = 0.1 case for both Markovian (Fig. 2) and non- 
Markovian case (Fig. 8) is compared, one finds that even 
at non-zero temperatures, non-Markovian baths works 
better in tripartite entanglement storage (Markovian case 
Af (t) dies at 1.25« -1 while non-Markovian baths sustains 
A f(t) at least till 2.5W” 1 while k = w” 1 ). Furthermore, 
we checked the same behavior extends down to S ^ 0 and 
lower oll values cases as well (not shown here). 
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(a) (b) 




(c) 



FIG. 8: Finite temperature effect on non-Markovian tripartite entanglement. Low temperature regime with 
N = 0.1 and the rest of parameters same as in Fig. 3 with time measured in w -1 units. As expected, inclusion 
of finite temperature in non-Markovian baths resulted in a relatively fast decay of A f(t) in all cases. We notive 
that in finite temperatures case n(t) attains maximum value rather quickly than the corresponding zero 
temperature case which results in an earlier entanglement decay. 


IV. CONCLUSIONS 
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zero and non-zero temperatures) with Lorentzian and 
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times. Either cavities start off in a maximally entan¬ 
gled W or GHZ state, we demonstrated that by prob¬ 
ing the non-Markovian bath parameters (detunning, cou¬ 
pling strength and temperature) temporal profile of en¬ 
tanglement can be controlled and made more robust 
against effects of decoherence in comparison to the cor¬ 
responding Markovian baths (with or without the possi¬ 
bility of photon hoping among cavities). In Lorentzian 
type of models we found far detunned double Lorentzian 
model with small coupling strengths, whereas in ohmic 
type of environments super-ohmic baths with small cut¬ 
off frequencies are excellent candidates that can maintain 
the three-cavity entanglement for at least thrice longer 
times (almost 10k -1 ) compared to all other Markovian 
and non-Markovian baths. 

Along with the calculation of non-Markovian tripartite 
entanglement with different environmental models, the 
main novelty of this work is the application of NMQJA 
to understand the dynamical features of J\f(t) under non¬ 
zero detunning. We also substantiated that according to 
NMQJA, time-dependent decay rates are sufficient for 
the understanding of the entanglement collapse and re¬ 
vival. For the cases when time-dependent decay rate re¬ 
mains positive, an average value of n(t) can be used to 
describe the faster/slower decline of negativity as a func¬ 
tion of time in different non-Markovian models. 
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Appendix A: Equations of motions of density matrix 
elements 


Note that the equations outlined below can de¬ 
scribe the Markovian baths (with and without pho¬ 
ton hoping) as well as the non-Markovian baths (with 
and without zero temperatures). By setting the 
Ji(ui) = Ki/2n and Ni(ui) = Ni below equations de¬ 
scribe the Markovian bath situation. While taking 
the time-dependent functions a*(f) = 0 and (3i(t) = 
fo dti / 0 °° Ji(w)e l f w-Wci ^ t-tl ^dw we enter into the non- 
Markovian bath regime with baths at zero tempera¬ 
tures. Photon hoping can be removed by considering 
£12 = £23 = 0. Upper dash represents time derivative of 
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matrix elements. 


pii(6 = 2Re(/3) (p 22 (f) + P3a(t) + p 44 (f)) - 6Re(a)pn(f) 

(Al) 

i°i 2 (*) = Pi 2 {t) (a — fi* — 6Re(a)) + 2Re(/3) (/ 935 (f) + /046(f)) 

(A2) 

/ 013 (f) = * 62 / 012 (f) + * 63 / 014 (f) + / 013 (f) (a - 6* - 6Re(a)) 
+2Re(/3) (p 2 s(f) + P 47 (f)) (A3) 

/oi 4 (f) = /°i 4 (f) (a - /3* - 6Re(a)) + 2Re(/3) (p 26 (f) + / 037 (f)) 

(A4) 

/ 015 (f) = / 015 (f) (2a - 2/1* - 6 Re(a)) + 2Re(/3)p 48 (f) 

(A5) 

pie (6 = * 63 /^ 15 (f) + * 62 / 017(6 + /oie(f) (2a - 2/3* - 6 Re(a)) 
+2Re(/3)p 38 (t) (A 6 ) 

/ 017 (f) = / 017 (f) (2a - 2/3* - 6Re(a)) + 2Re(/3)p 28 (f) 

(A7) 

Pis(f) = Pis(f) (~3a* - 3^*) (A 8 ) 

P 2 i (6 = P 2 i(f)(—a - /3 - 4Re(a)) + 2Re(/3) (// 53 (f) + Pe 4 (f)) 

(A9) 

P 22 (6 = *62 (p 2 3 (f) - / 032 (f)) + p 22 (f)(-4Re(a) - 2Re(/3)) 
+ 2 Re(a)/»n (f) + 2 Re(/ 3 )p 5 s (6 + 2R.e(/3)/»66(f) (AlO) 

P 23 (f) = *63P2 4 (f) 

-* 62 ^ 22(6 + P23(f)(-4Re(a) - 2Re(/3)) + 2 Re(/ 3 )/ 9 67 (f) 

(All) 

P 24 (6 = p 24 (f)(—4Re(a) - 2Re(/3)) + 2Re(/3)p 57 (f) 

(A12) 

P 25 (f) = / 025 (f)(a + /3 - 4R.e(a) - 4Re(/3))+ 

2 Re(a)/»i 3 (f) + 2Re(/3)p6 8 (f) (A13) 

P 26 (f) = * 63^25 (f) + * 62 ^ 27(6 — *6 2 P36(6 + P 26 (f)(a + /3 
—4Re(a) - 4Re(/3)) (A14) 

+2Re(a)pi 4 (f) + 2Re(/3)p 58 (f) 

/ 027 (f) = P 27 (f)(a + /3 — 4Re(a) — 4Re(/3)) (A15) 

P 28 (6 = /028(f) (2/3 - 2a* - 6Re(/3)) + 2Re(a)p 17 (t) 

(A16) 

P 31 (f) = P 3 i( 6 (-« - P - 4Re(a)) + 2Re(/3) (p 52 (f) + p 74 ( 6 ) 

(A17) 

P 32 (6 = P32(f)(-4Re(a) - 2Re(/3)) + 2Re(/3)p 76 (f) 

(A18) 

P3 3 (f) = *6 2 P23 (f) + *6 2 P32(6 + *63P3 4 (f) ~ *63P43(6 

+P 33 (f)(—4Re(a) - 2Re(/3)) + 2Re(a)pn(f) 

+2Re(/3)p 55 (6 + 2Re(/3)p 77 (f) (A19) 


P 3 4 (f) = P 3 4 ( 6 (—4Re(a) - 2Re(/3)) + 2Re(/3)p 56 (f) 

(A20) 

P 35 (f) = P 35 (f)(a + /3 - 4Re(a) - 4Re(/3)) + 2Re(a)pi 2 (f) 
+2Re(/3)p 78 (f) (A21) 

P36(f) = — *6 2 P26(f) + *63P35 (f) + *6 2 P37(f) + *63P46(6 + 

p 36 (f)(a + /3 — 4Re(a) — 4Re(/3)) (A22) 

// 37 (f) = p 37 (f)(a + /3 — 4Re(a) - 4Re(/3)) + 2Re(a)pi 4 (f) 
+2Re(/3)p 58 (f) (A23) 

P 38 (f) = P 3 s(f) (2/3 - 2a* - 6Re(/3)) + 2Re(a)p 16 (f) 

(A24) 

P 4 i(f) = P 4 i( 6 (—a - /3 - 4Re(a)) + 2Re(/3) (p 62 (f) + p 73 ( 6 ) 

(A25) 

P 42 (f) = P 42 (f)(—4Re(a) - 2Re(/3)) + 2Re(/3)p 75 (f) 

(A26) 

P 43 (6 = — *63P33(f) + * 6 2 p 42 (f) + *63P 44 (6 + 
p 43 (f)(—4Re(a) - 2Re(/3)) + 2Re(/3)p 65 (f) (A27) 

P 44 (f) = p 44 (f)(—4Re(a) - 2Re(/3)) + 2Re(a)pn(f)+ 


2Re(/3 )p 66 (f) + 2Re(/3)p 77 (f) (A28) 

P 45 (f) = P 45 (f)(a + /3 - 4Re(a) - 4Re(/3)) (A29) 

P46(f) = — *63P36(f) + *63P45(f) + *6 2 P47(6 
+P 46 (f)(a + /3 — 4Re(a) — 4Re(/3)) (A30) 

+2Re(a)pi 2 (f) + 2Re(/3)p 78 (f) 
p 47 (f) = P47(6(a + /3 - 4Re(a) - 4Re(/3)) 

+2Re(a)pi3(f) + 2Re(/3)p6 8 (f) (A31) 


p 48 (f) = P 4 s(f) (2/3 - 2a* - 6Re(/3)) + 2Re(a)p 15 (f) 

(A32) 

P 51 (*) = P 5 i(f)(- 2 a - 2/3 - 2Re(a)) + 2Re(/3)p 84 (f) 

(A33) 

P 5 2 (6 = P 52 (f)(—a - /3 - 2Re(a) - 2Re(^))+ 

2Re(a)p 3 i(f) + 2Re(/3)p 86 (f) (A34) 

P 53 (f) = — *63P36(f) + *62P52(6 + *63P54(6 
+P 5 3 ( 6 (—a - /3 - 2Re(a) - 2Re(/3)) (A35) 

+2Re(a)p 2 i(f) + 2Re(/3)p 87 (f) 
p ' 54 (6 = P 5 4 (f)(-a - p - 2Re(a) - 2Re(/3)) (A36) 

P 55 (f) = P55(6(-2Re(a) - 4Re(/3)) + 2Re(a)p 22 (f) + 

2Re(a)p 33 (6 + 2Re(/3)p 88 (t) (A 37 ) 

P56(f) = *63P55(f) + *62P57(f) - *63P66(6 + 
P56(6(-2Re(a) - 4Re(/3)) + 2Re(a)p 34 (t) (A38) 

P 5 7 (f) = P 57 ( 6 (~ 2 Re(a) - 4Re(/3)) + 2Re(a)p 24 (f) 

(A39) 

P 5 s(f) = P58(f) (/3 - a* - 6Re(/3)) + 2Re(a) (p 26 (f) + P 37 (f)) 

(A40) 

Pei (6 = Pei(f)(-2a - 2/3 - 2Re(a)) + 2Re(^)p 83 (f) 

(A41) 

Pe 2 (f) = P 62 (f)(—a - /3 - 2Re(a) - 2Re(^))+ 

2Re(a)p 4 i(f) + 2Re(/3)p 85 (f) (A42) 

P63 (f) = *63P53(f) + *62P62(f) + *63P6 4 (f) 
-*62P73(6 + P63(f)(-a - /3 - 2Re(a) - 2Re(/3)) 

(A43) 



12 


p'eiit) = P64(*)(-« - /3 - 2Re(a) - 2Re(/3))+ 

2 Re(o)/92i(*) + 2 R e{p)p 87 (t) (A44) 

P&S) = P65(i)(-2Re(a) - 4Re(/3)) + 2R e(a)p 43 (t) 

(A45) 

P(S6 W = — *^23P56(i) + i^23p65(t) + *£l2P67(^) — 

i£i2p76(t) + /0 6 6(i)(-2Re(a) - 4Re(/3)) (A46) 

+2Re(a)p 22 (f) + 2Re(a)p44(f) + 2Re(/3)p 88 (i) 

PctW = P67(i)(-2Re(a) - 4Re(/3)) + 2Re(a)p 23 (t) 

(A47) 

P6s(*) = P6s(t) (/3 - a* - 6Re(/3)) + 2Re(a) (p 25 (t) + P4 7 {t)) 

(A48) 

Pn(i) = P7i(i)(-2a - 2/3 - 2Re(a)) + 2Re(/3)p 82 (f) 

(A49) 

P 72 (^) = P 72 (t)(—a — /3 — 2Re(a) - 2Re(/3)) (A50) 

073 W = —i£l2P63(t) + i£l2P72{t) + *^23074(1) + 
/ 073 (*)(-a-/ 3 - 2 Re(a)- 2 Re(^)) (A51) 

+2Re(a)/9 4 i(i) + 2Re(/3)p 85 (t) 
p 74 {t) = p 74 (t)(—a — /3 - 2Re(a) - 2Re(/?))+ 

2Re(a)p 3 i (*) + 2Re(/3)p 86 (i) (A52) 

P 75 (*) = P 75 OX-2Re(a) - 4Re(/3)) + 2Re(a)p 42 (f) 

(A53) 

076^) = - *£12066^) + *£l2075 W + *£l2077(^) + 
P76(i)(-2Re(a) - 4Re(/3)) + 2Re(a)p 32 (f) (A54) 

077 (*) = p 77 (t)(-2Re{a) - 4Re(/3)) + 2Re(a)p 33 (f) + 
2Re(a)p44(t) + 2Re(/3)p 88 (t) (A55) 

P 78 (t) = P7s(t) (/3-a*- 6Re(/3)) + 2Re(a) (p 35 (t) + p 4 e(t)) 

(A56) 

Psi(t) = (-3a - 3(3)p 81 (t) (A57) 

082 W = 082 (t) (-2a - 2/3 - 2Re(/3)) + 2Re(a)p 7 i(f) 

(A58) 

083(*) = i£,i 2 P% 2 {t) + i&3P8i(t) + p 83 (t)(-2a - 2/3 - 2Re(/3)) 
+2Re(a)p 6 i(i) (A59) 

084 (*) = Ps4{t)(-2a - 2/3 - 2Re(/3)) + 2Re(a)p 5 i(t) 

(A60) 

P 85 (i) = P85(*)(-a - /3 - 4Re(/3)) + 2Re(a) (p 62 (t) + P 73 W) 

(A61) 

PseW = i& 3 p 85 (t) + i£i 2 p87(0 + pm(t)(-a - /3 - 4Re(/3)) 
+2Re(a) ( 0 . 52 (f) + p 7 i(t)) (A62) 

0s 7 ( i ) = 087(f) (-a - (3 - 4Re(/3))+ 

2Re(a) ( 053 (f) + 064 (f)) (A63) 

P 88 (*) = 2Re(a) (p 5 5 (*) + 066(f) + 077 (f)) - 6Re(/3)p 88 (f) 

(A64) 


function method introduced by Piilo et. al. In stan¬ 
dard Markovian QJA we think of actual/fictitious photo 
detectors surrounding the system of interest such that 
whenever a photon is lost by the system any one of the 
detectors register it by making a click. According to the 
QJA then in any given small time interval we have one 
of two situations: 

(I) Jump Case: A quantum jump takes place, and 
we apply the output (annihilation) operator Jj ([30]) as¬ 
sociated with the jth jump channel/detector (which in 
this context is also termed the “jump operator”) to the 
state of the system. 

(II) No Jump Case: No detector clicks, the sys¬ 
tem evolves according to the following non-Unitary 
Schrodinger equation: 


ih*—— 

dt 


4f(f)) = Hnh 4'(f) 


(Bl) 


where the “Hamiltonian” appearing in this equation is a 
non-Hermitian operator, which is the sum of two parts. 
The first part is Hermitian and is given by the system’s 
Hamiltonian, the second part is anti-Hermitian and is 
constructed from the jump operators. In total we have 

Hnh = H sys — — ^^ 7 jJjJj (B2) 

3 

7 j is the rate at which excitation decays through the jth 
channel. Note that for the sake of clarity here instead of 
using k we are using 7 as decay functions. Based on the 
cycle of these oscillations jump part of the non-Hermitian 
Hamiltonian is then divided into negative and positive 
parts as: 




3 + 


3~ 


(B3) 

Note that corresponding to positive and negative decay 
rates positive and negative jump operators Jj± 

are introduced. Interpretation of these negative and pos¬ 
itive jump channels is as follows: 

During the positive channel/cycle of 7 j(t) ( 7 j(t) > 0), 
system behaves as the source of the excitation and the 
environment becomes the target of the information flow. 
Hence after the positive jump system loses information 
thus producing the state reduction in an abrupt way as: 


I’M 


I’M = 


Jj+ |4 f a) 

I MM I 


(B4) 


Appendix B: Review of non-Markovian Quantum 
Jump Approach 

Closely following the reference [52], in this section we 
briefly review the NMQJA/quantum Monte Carlo wave 


subscript a is introduced to mention that at time t there 
are N a number of ensemble members in the source state 
|4' a (f)}. |'I' Q ,/(t)) is the target state in this scenario with 
N a i members. Occurrence of jump is a random process 
which in the positive decay cycle happens in a small time 
interval [f, t + St] with the probability: 

P?{t) = J j+ (t)(* a (t)\j} + j j+ \* a (t))6t (B5) 
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Note that the Markovian QJA is a special case of this pos¬ 
itive decay cycle situation with positive time-independent 
decay rates. 

During the negative cycle/channel ( 7 j-(t) < 0) jump 
process reverses it’s direction i.e. memory-full environ¬ 
ment becomes the source and system acts as the target 
of the information 

I’M <— I’M = ( B6 ) 

\\Jj- I’M II 

It is worth mentioning here that click description of 
quantum jumps is (strictly speaking) valid for positive 
jumps only. For negative jumps, the click picture is 
still an under debate issue (for further details see for in¬ 
stance [42, 43]). During negative cycle, probability of 
negative/non-Markovian jump turns out to be: 

MM) = ^|7M)KM*)|Jj_ Jj-lVa'iW 

(B7) 


Note that the negative jumps are actually responsible 
for producing memory effects (non-Markovianity) in the 
system dynamics. In above equation negative jump prob¬ 
ability is weighted by the ratio of number of ensemble 
members in the target/system state to the number of 
members in the source/environment state. This implies 
that negative jump probability is going to be zero if there 
are no members in the target state. Moreover we have 
taken the absolute value of negative decay rate to ensure 
that the probability remains non-negative. Sign of the 
negative decay (as already mentioned) describes the re¬ 
verse flow of information and therefore non-Markovian 
quantum jump approach doesn’t suffer the issue of nega¬ 
tive probabilities [28]. Finally note that one can take the 
ensemble average over all different possible trajectories 
(no, positive and negative jump situations) to construct 
the full density matrix which will be equivalent to the 
density matrix obtained by directly solving the full non- 
Markovian Master equations [22] , 
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